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Abstract 

The Nikolaevskiy equation was originally proposed as a model for seismic waves and is also a 
model for a wide variety of systems incorporating a neutral "Goldstone" mode, including electro- 
convection and reaction-diffusion systems. It is known to exhibit chaotic dynamics at the onset of 
pattern formation, at least when the dispersive terms in the equation are suppressed, as is commonly 
the practice in previous analyses. In this paper, the effects of reinstating the dispersive terms are 
examined. It is shown that such terms can stabilise some of the spatially periodic traveling waves; 
this allows us to study the loss of stability and transition to chaos of the waves. The secondary 
stability diagram ("Busse balloon") for the traveling waves can be remarkably complicated. 
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I. INTRODUCTION 

In 1989, Nikolaevskiy [1] derived a model for longitudinal seismic waves, in the form 
of a one-dimensional partial differential equation for a displacement velocity. Although 
Nikolaevskiy' s equation included dispersive terms, most subsequent analysis has treated a 
simplified version of the PDE, in which these terms are omitted. This reduced form is now 
generally known as the Nikolaevskiy equation, which may be written in the form 
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where r is a control parameter. The equation ([I]) has been proposed as a model for several 
other physical systems, including phase instabilities in reaction-diffusion equations j^j , elec- 
troconvection [3J and transverse instabilities of fronts ^]. More generally, (JTJ can be regarded 
as a simple model of a pattern-forming system with an instability at finite wavenumber and 
a neutral "Goldstone" mode arising from symmetry {3, 5]. 

The uniform state u = of (JTJ becomes unstable at r = to spatially periodic "roll" 
solutions, with wavenumbers around k = 1. However, these, in turn, are themselves all 
unstable at onset in sufficiently large domains {3]; this unusual instability arises from the 
neutral mode at wavenumber k = 0. In fact, numerical simulations show that the Niko- 
laevskiy equation exhibits spatiotemporal chaos at onset a,l6|. The scalings associated with 

of the equation has stimulated significant investigation [8J. 

Although in some applications (such as the instability of fronts (4j) the omission of dis- 
persive terms is justified on symmetry grounds, this is not the case in the original context 
of a model for seismic waves 

Earlier work that has considered the effects of dispersion includes the paper of Mal- 
omed 9)], who reinstated one dispersive term in the Nikolaevskiy equation and analysed 
the secondary stability of traveling-wave solutions by means of coupled Ginzburg-Landau- 
type equations for the amplitude of the traveling waves and a large-scale mode. His results 
showed that dispersion could stabilize waves; however, his derivation was not entirely asymp- 



totically self-consistent [3[. Kudryashov and Migita [10] showed, on the basis of numerical 
simulations, that traveling waves can be stabilized by the presence of dispersive terms in the 
Nikolaevskiy equation. It is also known that in the related Kuramoto-Sivashinsky equation, 
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the introduction of a dispersive term can stabilize periodic traveling waves 11]. 

Our aim in this paper is to provide a systematic examination of the effects of dispersion. 
By varying the parameters corresponding to dispersion, we can find when dispersion stabi- 
lizes traveling waves and investigate how the chaotic state in the non-dispersive equation 
arises as the dispersion is reduced. 

In the following section we give the form of the equation and the traveling waves under 
consideration. Computational results on the stability of these waves are given in Sec. III. 
The stability analysis of the waves is complicated and depends on the magnitude of the 
dispersion terms; three different scalings are considered in Sees. IV, V and VI. Sec. VII 
illustrates some numerical simulations of the Nikolaevskiy equation with dispersion, and our 
conclusions are summarized in Sec. VIII. 



II. THE NIKOLAEVSKIY EQUATION WITH DISPERSION 



We examine the Nikolaevskiy equation with dispersion in the form 
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where a and /3 are the dispersion coefficients. This equation is thus the one originally 
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proposed by Nikolaevskiy [l| (and later examined in |9|, [10|), with all spatial derivatives up 
to the sixth appearing on the right-hand side. In the numerical simulations presented in 
Sec. I VI I| we shall impose the periodic boundary condition 



u(x + D,t) = u(x, t) 



(3) 



for some domain length D. 

Before proceeding, we note that fl2]) has the same Galilean symmetry (x i— > x + Vt, 
u h> u + V) as the nondispersive equation ([1]). Moreover, in view of the Galilean symmetry 
and the observation that, when the boundary condition is imposed, 

d r D 

u(x, t) dx = 0, 



dt J0 

the spatial average of u may be set as zero (by transforming to a moving frame of reference 
if necessary). The reflection symmetry (x h-> —x, u \- > —u) of (CQ) is broken by the presence 
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FIG. 1: Plot of A r for the case r = 0.1: note the linearly growing modes with wavenumber around 
k c = 1, and the weakly damped large-scale modes close to k = 0. 

of the dispersive terms. However, there is a symmetry x \- > —x, u \- > —u, a h-> —a, (3 h-> — /3; 
as a consequence of this symmetry we need consider only the case > 0. 
Linearization around the steady state w = yields the dispersion relation 

A = k 2 [r - (k 2 - l) 2 ] + ik 3 (k 2 (3 - a) 

for Fourier modes proportional to e lhx+xt . Thus in general these perturbations take the form 
of traveling waves, with phase speed 

c P = = k 2 (a - k 2 (5) (4) 



and group velocity 



Cg = -^ = k 2 (3a-5f3k 2 ). (5) 



The real part of the growth rate, A r , is plotted in Fig. [fl for r just above the threshold value 
r c = for the onset of instability. This figure shows that there exists a band around the 
critical wavenumber k c = 1 of linearly growing modes, and a neutral mode at k = (the 
so-called "Goldstone mode"), which significantly affects the nonlinear dynamics of (J2]). 

Just beyond the onset of instability of the zero solution, it is straightforward to carry out 
a weakly nonlinear analysis of (j2J), with 

r = e 2 r 2 . (6) 
This analysis reveals that there are traveling- wave solutions of the form 

u ~ eaoe^*-^ + c.c, (7) 



where the wavenumber k = 1 + eq. The amplitude turns out to be given by 

a = 6(r 2 - 4g 2 ) 1/2 (l + U a ~ ^ff 1 " ( 8 ) 
and the speed of the wave is 

s = c p - ±e 2 (r 2 - V)(a - 5/3) + o(e 2 ), (9) 

where c p is given by (jl]) and the second contribution to s reflects (weakly) nonlinear effects. 
So, regardless of the values of the dispersion parameters a and (3, such spatially periodic 
solutions exist for r 2 > 4g 2 . We now turn to the question of the secondary stability of these 
solutions. 

III. SECONDARY STABILITY OF TRAVELING WAVES: NUMERICAL RE- 
SULTS 

In this section we first outline a numerical method for the calculation of the nonlinear 
traveling waves and their secondary stability, and then give the results of these computations, 
showing the stability boundaries of traveling waves. 

A. Numerical method for calculating secondary stability 

To calculate the secondary stability of a traveling wave solution for given values of the 
parameters, we first find the traveling wave solution u(x, t) = f(z), where z = x — ct. Here, 
c is the nonlinear wave speed, which in general is not exactly equal to the linear wave speed 
c p 01]). We approximate the solution numerically using the truncated Fourier series 

N/2 

f( Z )= Yl u ^ nkz . 

-N/2+1 

Substitution in fl2]) (and calculation of the nonlinear term pseudospectrally) yields a system 
of nonlinear equations (solved in Matlab) for the Fourier coefficients of f{z), together with 
c, which is determined from 

c ( D {f) 2 dz = a / V) 2 dz-fi [ D (f") 2 dz + [ D f{f) 2 dz, (10) 
Jo Jo Jo Jo 
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where D = 27r/k is the length of the domain and k is the wavenumber of the solution under 
consideration. The expression f flQj) follows from multiplying (J2J) by f'(z) and integrating over 
the domain, using integration by parts multiple times. To compensate for the additional 
unknown c, we have an additional equation from the fact that we may choose the phase of 
the wave, for example by specifying that U\ is real. 

After calculating the solution, we construct the eigenvalue problem for perturbations. 
If we suppose that u(x,t) = f(z) + u(x,t), then substitution in (T5]) yields the linearized 
perturbation equation 



du d 2 



dt dx 2 

We take 



1 + 



d 2 x 2 



dx' 
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N/2 



u = e at+ipz v ne inkz , 



-N/2+1 

where all possible eigenfunctions may be captured by limiting consideration to — k/2 < p < 
k/2. The resulting eigenvalue equations to determine the growth rate a are then 

N/2 N/2 
-N/2+1 -N/2+1 

where C = K 2 (r — (1 — K 2 ) 2 ) — iaK^ + if3K^ and K n = p + nk. The eigenvalues of this 
system are computed numerically. By examining the largest real part of all eigenvalues a 
for a large sample of values of p in the relevant interval, we determine whether the original 
traveling waves are stable or unstable. In the following section we provide some stability 
diagrams based on the above method. 

In determining our results, we have been careful to check that: adequate samples in p 
are taken (too few, particularly for small values of p, can lead one to miss certain small 
regions of instability); adequate Fourier modes are taken in determining both the original 
solution and the perturbations; adequate samples are taken in parameter space to determine 
all regions of stable rolls. Typically, 300 values of p are used, with N = 16. 



B. Results 



Now we present the secondary stability diagrams. The first case considered here is setting 
(5 = and varying a — see Fig. [2j When a is small (a = 1/2), there is a very small region 
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FIG. 2: The secondary stability regions of traveling waves of ([2]), calculated numerically for (a) a = 
1/2, (b) a = 2 and (c) a = 5, all for (3 = 0. Shown are the marginal curve r = (1 — k 2 ) 2 (solid line) 
and the secondary stability boundary of the traveling waves (dashed line), with stability between 
the dashed lines. 
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FIG. 3: The secondary stability regions of traveling waves of ([2]), calculated numerically for (a) f3 = 
5, (b) j3 = 5.5, both with a = 40. Shown are the marginal curve r = (1 — k 2 ) 2 (solid line) and 
the secondary stability boundaries of the traveling waves (dashed line), with stability between 
the dashed lines. To clarify the regions of stability/instability, the "s" indicates one of the stable 
regions. 

of stable waves in the (k, r) plane. The stable region is a thin strip, confined to small values 
of r; in this case, for r > 0.0078 all rolls are unstable. 

For larger a this strip of stable waves is longer and wider; for example at a = 2 there are 
some stable rolls up to r « 0.22, and at a = 5 the stability region extends at least as far 
as r = 0.9. Furthermore, it is apparent for a = 5 that a symmetrical Eckhaus-like stability 
region is present for very small values of r (from the numerical results themselves, it seems 
to be present in all three cases, but is visible only in the last plot of Fig. [2]). The shrinkage 
of the region of stable traveling waves for small a is consistent with there being no stable 
rolls at all in the nondispersive case. 

While an exhaustive examination of the secondary stability diagrams across (a, 0) pa- 
rameter space is infeasible, it is worthy of note that these diagrams may be extremely 
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FIG. 4: The secondary stability of traveling waves of ([2]) calculated numerically for (a) fixed (3 = 
and r = 0.01 in (k,a) parameter space, (b) fixed a = 40 and r = 0.1 in (k,/3) parameter space. 
The marginal curve is represented by the solid lines; traveling waves are stable inside the dashed 
lines. 

complicated. A good example arises if we set a = 40 and vary /3 — see Fig. [3J For (3 = 5 
there is a small Eckhaus-like stability region for r < 0.001; for larger values of r, there re- 
mains a single stability region. For larger /3, however, the stability region splits into several 
parts; for example, at /3 = 5.5 there may be up to five separate intervals of stable traveling 
waves for a given value of r. 

Above we have presented our secondary stability diagrams in the (k, r) plane, for fixed 
values of a and (3. If our interest is in the effects of dispersion on the stability of traveling 
waves then it is more instructive instead to fix r and present results in either the (k, a) or 
the (k,/3) plane. Our first example is for r = 0.01 and (3 = — see Fig. H^a). Given this 
value of r, the traveling waves exist for 0.9487 < k < 1.0488. We expect that if a is small 
enough then all roll solutions are unstable; this is indeed the case. For larger values of a, a 
region of stable rolls appears. In Fig. 11(b), we show a second case, where we fix a = 40 and 
r = 0.1, to emphasize that the structure of the stability region may be rather complicated, 
exhibiting a sensitive parameter dependence. 

IV. SECONDARY STABILITY OF TRAVELING WAVES: a, (3 = 0(1) 

In this and the following two sections, we analyse the secondary stability of traveling waves 
(J7J). The most straightforward case arises when the dispersion parameters a and (3 are each 
0(1). To contrast with later sections, we shall characterize this case as strong dispersion. 
Whereas the nondispersive Nikolaevskiy equation has no stable spatially periodic states, 
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Kudryashov and Migita [10(| found stable periodic waves in their numerical simulations of 
the dispersive PDE ([2]), in this regime. 

We begin by introducing the weakly nonlinear expansion 

u = eui + e 2 u 2 + e 3 w 3 H , (12) 

with r given by (J6]). Then substitution in ([2]) and consideration of successive orders in e 
leads to the following. 
At 0(e), we find that 

Ul = Ae i{x ~ Cot) +c.c, 

where Co = a — (3, and where the amplitude A varies slowly in space and in time, in principle 
depending on the slow variables 

X = ex, t = et, T = e 2 t. 

A consideration of the terms proportional to e 1 ^ 00 *' at 0(e 2 ) then shows that in fact A = 
A(£, T), where 

£ = X - (3a - 5/3)t = X - vr 

is a coordinate moving at the group velocity of the waves. Then solving the problem at this 
order in e yields 

"^- 36( 1+1 (f-5, fl )/6) e ^" + - + / - 
Here / is a slow varying function of X, r and T, chosen to appear at this order to balance 
forcing terms appearing at the next order in e. 

At 0(e 3 ), we find, from the respective consideration of the terms in ((2]) proportional to 
e i(x-c t) an( j e oi(x-cot)^ amplitude equations 

I - 

Since A = A(£,T), the second amplitude equation suggests taking / = /(£, T), in which 



case p4j) becomes 

9/ _ 9|A| S 



— 1> 



<9£ ^ 
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FIG. 5: Diagram showing the sign of 1 + a(b + d) in a/3 parameter space. Regions with "s" are 
where 1 + a(b + <£) > 0, so that a limited band of plane waves is stable, as in (JEJ); a "u" indicates 
where 1 + a(b + d) < 0, and all plane waves are unstable. 

and hence vf = \A\ 2 + K(T), for some K(T). However, the constraint that the spatial 
average of u should be zero gives K(T) = — (\A\ 2 ), where the angle brackets denote the 
average in £. Thus 

f= dm±w (15) 

and the amplitude equation (fT3l) becomes the nonlocal Ginzburg-Landau equation 

9A ( 1 -i(«-5/3)/6, „ l2 (\A\ 2 )-\A\ 2 \ A ft , d 2 A , . 



It is worth mentioning that in view of ([15]) the present scaling breaks down when v is small; 
in particular, this is the case when a and (3 are both small, and this case will be considered 
in later sections. 

It is helpful in analysing ( JT6l) to put it in canonical form by rescaling all the variables, to 
give 

8 A f) 2 A 

— = A + id((\A\ 2 ) - \A\ 2 )A + (1 + ia)— - (1 + ib)\A\ 2 A, (17) 



where 



3a -10/3 , 5/3 -a , 36 + (5/3 - a) 2 
a = ■ , o = — - — , a - 



4 6 ' v 

Equations similar to (fl?)) . including a nonlocal nonlinear term have been derived and stud- 



ied in the context of convection in a rotating annulus 12[ and in electrical and magnetic 
systems HQ. 
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Equipped with (|T7|) . we are now in a position to explore the secondary stability of weakly 
nonlinear spatially periodic solutions of the dispersive Nikolaevskiy equation. Such solutions 
correspond to plane- wave solutions of (fiT|) . which exist in the form A = Pe^ ulT+q ^\ with 
P = (1 — g 2 ) 1 / 2 and to = q 2 (b — a) — 6. To study the stability of the plane- wave solution, we 
write A = (1 + p(£, T))Pe l ^ T+q ^\ which, after substitution in ( ITTj) and linearization in the 
perturbation p, yields 



^ = (1 + ia) + 2ig^ J - (1 + ib)P 2 (p* + p) + idP 2 ((p + p*) - (p + p*)) . 

Then upon setting p(£, T) = R(T)e lL ^ + S*(T)e _li ^ and equating the coefficients of e lL ^ and 
e~ lL ^, we have 

^ = -(1 + ia)L(LR + 2qR) - (1 + i6)g 2 (i? + S) - idQ 2 (R + 5), 

A Q 

— = -(1 - ia)L(L,S - 2g5) - (1 - i6)g 2 (i? + 5) + icZg 2 (i? + S). 

Finally, with R(T) and S(T) proportional to e^ T , and expanding the growth rate in powers 
of the perturbation wavenumber L, we have the dispersion relation 

fj, = -2iq(a-b-d)L + L 2 p- 2 (-1 - a(b + d) + q 2 [3 + 2(6 + df + a(6 + d)])+0(L 3 ). (18) 

If we suppose (as is generally the case) that a ^ b + d then it is apparent from ffl8|) that 
the solution has a long-wavelength oscillatory instability whenever 

q 2 (3 + 2(6 + d) 2 + a(6 + d)) > 1 + a(6 + d). 

Since 1 + a (6 + <f) < 3 + a (6 + d) + 2(6 + d) 2 , we see that stability is determined by the 
following. If 1 + a(b + d) > then 3 + a(b + d) + 2(6 + d) 2 > 0, and the plane- wave solutions 
are stable provided 

9 9 1 + a (b + d) _ . 

^" 2 <"^ WTWTWT^+d) <h (19) 

If instead 1 + a(b + d) < 0, then all plane waves are unstable. We note that setting a = 6 = 
d = reduces (TTT|) to a real Ginzburg-Landau equation for A, and our results reduce to the 



usual Eckhaus instability (with stability for q 2 < 1/3) 15l. Il6|. 

To apply this result to ([2]) it is necessary to indicate the regions in a, ft parameter space 
in which the quantity 1 + a (6 + d) is positive or negative. In Fig. [5j regions denoted by "s" 
indicate where 1 + a(b + d) > 0, so that there are some stable plane waves, as in ffl9l) ; those 
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FIG. 6: Numerical simulations of the amplitude equation (|17|) : in each case the real part of A is 
plotted as a function of £ and T. In (a) and (b) a = 10, whereas in (c) a = 8.4; in each case (3 = 2.6. 
The initial condition in each case is a plane wave, with n wavelengths in the computational box 
— 32-7T < £ < 32-7T, plus small-amplitude random noise: (a) n = 28 (hence q = 0.875); (b) n = 10 
(q = 0.3125); (c) n = 20 (q = 0.625). 

regions denoted by "u" show where 1 + a(b + d) < 0, and hence all plane waves are unstable. 
As discussed in Sec. HI1 only the region (3 > need be presented. The existence of a stable 
region in Fig. |5]is consistent with the numerical results of Sec. IIII( for example Fig. [2] shows 
a stable region when a = 0(1) and (3 = 0. 

In Fig. [HI we illustrate the considerations above with some numerical simulations of the 
modified complex Ginzburg-Landau equation ( ITT)) . Our numerical code is pseudospectral, 
and uses exponential time differencing jl7| . In each case the initial condition is a plane wave 
plus small-amplitude random noise. For the simulations illustrated in Fig. E(a) and (b), 
1 + a(b + d) > 0. The two plots show the fate of initial conditions in the unstable and 
stable regions of Fig. El respectively. In each case, a stable plane wave is obtained at large 
T. Figure El^c) shows the development of instability in the case 1 + a(b + d) < 0, where all 
plane waves are unstable. Here the solution is persistently time-dependent. 

The analysis above tells us about the secondary stability of traveling-wave solutions of 
the dispersive Nikolaevskiy equation when a, (3 = 0(1), and the results are summarized in 
Fig. [5j We may think of this analysis as holding for any fixed a and (3 (not both zero) in the 
limit as r — > 0; thus we expect the lowest part of the secondary stability diagram in (r, k) 
parameter space to reflect Fig. [5j 

However, as indicated earlier, when a and (3 are both small, the analysis above does 
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not hold, and requires reconsideration. We should expect such analysis to break down 
in this limit, because Fig. [5] is inconsistent with the known behavior of the nondispersive 
Nikolaevskiy equation (a = (3 = 0), for which all rolls are unstable at onset {3, 5|. Thus in 
the next section we consider smaller values of a, (3. 



V. SECONDARY STABILITY OF TRAVELING WAVES: a,P = 0(e 3 / 4 ) 

It turns out, after some experimentation, that small a and first lead to a new scaling if 
we adapt the scaling first used by Tribelsky and Velarde [3[ for the nondispersive case, and 
extended by Cox and Matthews jj] to a damped version of the Nikolaevskiy equation. In 
this scaling the original traveling waves remain 0(e), but the perturbation to the traveling- 
wave amplitude is 0(e 3//2 ) and the large-scale mode is 0(e 7//4 ); furthermore, slow spatial and 
temporal variations of perturbations take place on scales given by X = e 3//4 x, T = e 3 / 2 £, 
r = e 3 / 4 t. (Note that these slow variables are different from those of the previous section, 
but our notation for slow variables is consistent within sections.) To allow the development 
of consistent amplitude equations for the perturbation we then take 

a = e 3 / 4 d, (3 = e 3 / 4 /3. 

Applying a weakly nonlinear analysis to (T5]) gives 

u = e (a + e 1 / 2 a(X, T))e iM + c.c. + e 7 / 4 /(X, T) + • • • , (20) 

where ao = 6^r 2 — 4g 2 , 

M = (1 + eq)x -CT- e 1/4 vqT + e l/A ifj(X, T), 

c = a — j3 and v = 3a — 5/3. Here a(X, T) represents disturbances to the amplitude of 
the pattern, ip(X,T) represents corresponding disturbances to the phase of the pattern 
and f(X,T) is a large-scale mode. Substitution of u, as given by ( )20l) . in fl2]) requires the 
consideration of the problem at successive orders in e 1 / 4 . After much consequent algebra, 
we find the (nonlinear) amplitude equations 

dip d 2 i/j A dip 

&T = dx 2 ~ f ~ v ~dx' 
fff ~ ox 2 a °~dx' 

da d 2 a ( dip\ 2 dip da 

&T = A dx 2 ~ Aao [dx ) - 8aoq dx~ v dx- 
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Note that dispersion is represented in these equations only through the terms vipx and vax, 
representing advection of the pattern envelope with the group velocity v. Note also that the 
group velocity of the large scale mode / is zero, and hence no corresponding term appears 
in the second of these equations. 

The three amplitude equations may be reduced to the single (nonlinear) phase equation 

d d 2 A d \ V d d 2 \ , 2 / dip \9 2 ^ 



#r ax 2 ox J \ot dx 2 ) r u v^x 7 ax 2 

Then linearising this equation and setting ip = e lLX + aT yields the dispersion relation 

a 3 + 9a 2 L 2 + 2AaL 4 - v 2 aL 2 + 16L 6 - v 2 L 4 - lQa 2 qL 2 + iv(2a 2 L + 10aL 3 + 8L 5 ) = 0. (22) 

Before considering this dispersion relation for general L, it is helpful to consider the two 
limiting cases, of small and large L. First, if L is small, then <r 3 ~ Wa^qL 2 . Thus, to leading 
order in L, a = (T2/3L 2 / 3 , where cr 3 ^ = IQa^q] hence all traveling waves are unstable if L i 



is 



small. On the other hand, if L is large, then we have a 3 + 9a 2 L 2 + 24aL 4 + 16L 6 w 0, and 
so cr ?y —L 2 or — 4L 2 (twice); hence traveling waves are stable to large-L disturbances. In 
summary, all traveling waves are unstable at onset (provided a^q 7^ 0; in fact we shall see 
later that when a^q is suitably small, we shall need to reconsider this conclusion). The rest 
of the section provides more details of the instability, for general values of L. 

In order to find the secondary stability boundary for the traveling waves, we set a = iQ 
in the dispersion relation (|22|) . where Q is real. From the real and the imaginary parts, we 
obtain 

^ 2 16 . 16 2 v 2 t2 10 A r _ 

Q 2 L H alq H L H vLQ = 0, 

9 9 oy 9 9 

tt 3 - 24ttL 4 + v 2 VLL 2 + 2vLtt 2 - 8vL 5 = 0, 

and then after eliminating Q between these two equations we find that this stability boundary 
is given by 

16a^g 3 - 2500L 12 + 2100L 8 a 2 ) g + 384L 4 a 4 g 2 - 200v 2 L 10 - Av 4 L 8 - Uv 2 L 6 a 2 q + v 2 L 2 a 4 q 2 = 0. 

(23) 

We note that in this equation L and v appear only as even powers and thus we can restrict 
our attention to positive L and v with no loss of generality. However, both even and odd 
powers of q occur, so no such economy is possible in considering q (indeed, in the light of j^|, 
we should expect different behaviors for q > and q < 0). 
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FIG. 7: Predicted secondary stability boundaries of spatially periodic solutions of the Nikolaevskiy 
equation, (a) Nondispersive case (v' = 0), (b) dispersive case (V = 5). 

For the case of no dispersion, Tribelsky and Velarde [3] showed that, according to the 
present scaling, there is monotonic instability of the rolls with q > (with unstable distur- 
bances having < L < (a^g) 1 / 4 ). By contrast, oscillatory instability occurs for rolls with 
q < (unstable modes having < L < (— 2agg/25) 1 / 4 ). 

It is convenient to present our results for the dispersive case in terms of the rescaled 
variables q' = q/r^ 2 , V = L/r^ and v' = v/r^ 8 . Figure [7] illustrates the regions of 
stability and instability of the traveling waves, in the cases v' — and v' = 5. Note that in 
the dispersive case all instabilities are oscillatory. 

We should view with caution the conclusion above that all traveling waves are unstable, 
because it relies crucially on the assumption that a^q is not small. The stability analysis 
above breaks down if q or Oq are small; the true stability properties of corresponding traveling 
waves will be investigated in the next section. 



VI. SECONDARY STABILITY OF TRAVELING WAVES: a,/3 = O(e) 

In this section, we investigate the cases of traveling waves with wavenumber close to k — 1 
or close to the marginal stability boundary, in other words those for which in the previous 
scaling a^q <^ 1. 



A. Traveling waves with close-to-critical wavenumber 

In order to resolve the secondary stability problem for traveling waves with wavenumber 
close to k c — 1, we set k — 1 + e 2 q, as was done for the nondispersive case by Tribelsky and 
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Velarde Then a distinguished balance occurs for a, f3 = 0(e); so we write 

a = eat, ft = eft. 

Upon setting c— a — 8, r = e 2 r 2 , v = 3a — 5ft, X = ex, r = et and T = e 2 t (the scalings for 
X and T being as in |3j), we find from (J2]) that 

u = e(6v^ + e 2 a(X, T))e iA/ + c.c. + e 3 f(X, T) + ■ ■ ■ , 

where now 

M = (1 + e 2 g)x -ct + e(-vq + |r 2 (d - 5ft))T + e^(X, T). 

The terms in M involving a and ft correspond to nonlinear effects of the finite traveling-wave 
amplitude on the speed of the waves; see (Q. 

After much algebra, the relevant (nonlinear) amplitude equations are found to be, at 
0(e 4 ) and 0(e 5 ), 

M- 4 ^i-f-^ (24) 
dT ~ dX* 1 V dX' [ ] 

dj_ _ i2r 1/2 — (25) 

dT ~ dX 2 Ur2 dX> {2b} 

da A d 2 a n A 1/a / dtp \ 2 A da i/ 2 9/ 



+ 6 ^(_S 5 + f r2 + 1 2^ + (^-3 & )^)§. ,26) 

We note that in these equations the influence of dispersion arises not only through the terms 
involving the group velocity v, but also through the term 10/3 — 3d in the equation for a?, 
in contrast to the previous case. 

To determine the stability of the traveling waves, these equations are linearized; for 
solutions proportional to e lLX+ °" T , we find the dispersion relation 

a 3 + 9cr 2 L 2 + 24aL 4 + 16L 6 + 528r 2 L 2 
- 576r 2 qL 2 + 82r 2 oL 2 - 568r 2 L 4 + 2r 2 cx 2 - v 2 oL 2 - LV 
+ i(2r 2 vaL + 360r 2 /3X 3 + 8vL 5 + WvaL 3 + 2va 2 L + 2r 2 vl?) = 0. (27) 

As in the previous section, in the limit of large L, all eigenvalues have negative real part. 
By contrast, in the limit of small L, if we expand a = 0\L + o 2 L 2 + • • • , then from (1271) we 
find that o\ satisfies 

r 2 o\ - 288r 2 q + 2<oAr\ + %r 2 vo x = 0, 
16 
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FIG. 8: The stability boundaries for traveling waves according to the amplitude equations (|24 
PB|). in the special case /?' = 0, for different values for a': (a) a' = 1, (b) a' = 3, (c) a' = 4, (d) 
a' = 5, (e) a' = 5.7 and (f) a' = 7. See the text for more details. 

whereas cr 2 is determined from 

a\ + 82r 2 ai + 4r 2 <Tior 2 - uVj + 2i(r 2 w 2 + ua? + r 2 v + 180r 2 /3) = 0. 
The first of these gives 

<ri = - (-iv ± v/-{) 2 + 1152g- 1056r 2 ) , (28) 

and so traveling waves are certainly unstable if their wavenumber satisfies q > llr 2 /12 + 
■0 2 /1152. The term -0 2 /1152 indicates that these waves become more stable with respect to 
this instability in the presence of dispersion. If instead q < llr 2 /12 + v 2 /1152, then <j\ is 
purely imaginary, and stability is determined by 



-72vq/r 2 + 171-0/2 - 180/3 91 72 
°"2 = ± 7= 1" 9> 



(29) 



^v 2 - 1152g + 1056r 2 2 r 2 
a consideration of which shows that these waves are made more unstable to the long- 
wavelength oscillatory instability in the presence of dispersion. 

Analysis of the stability boundaries to disturbances of general L is rather involved, and 
we do not present the details here. Furthermore, the parameter space is large enough to 
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preclude our making general statements; instead we consider some illustrative special cases. 

1 li 

To present the conclusions most generally, it is helpful to introduce q' = qjr^-, L' = L/r 2 , 
a' = a/r^ 2 , f3' = (3/r^ 2 and v' = vjr^ 2 . 

Let us begin by considering the special case /3' = 0. Figure [S] shows where traveling 
waves with different values of q' are stable and unstable to perturbations with wavenumbers 
L'\ each panel in the figure corresponds to a different choice of a'. In understanding the 
sequence of transitions in the topology of the various panels, it is helpful to first consider the 
behavior of the stability boundaries for small a' (and hence small v'), in particular in the 
U — limit. We have seen above that the right-hand stability curve (labeled R) intersects 
the q' axis at q' + = 11/12 + v' 2 / 1152. For small v, it follows from (|29p that the left-hand 
stability curve (labeled T) intersects the q' axis at q'_ ~ 91/144 + 5|t> 'I/26568 1 / 2 . Thus as 
a' is increased from zero, q'_ moves to the right more rapidly than does q' + . Eventually, 
at some sufficiently large value of a', q'_ = q' + , and all traveling waves are unstable in the 
limit V = 0. On the other hand, when a' is large, q'_ halts at q'__ = 131/144. However, q' + 
continues to increase, and this results in the appearance of a small-L' stability region. In 
fact, for sufficiently large a', some rolls are stable to disturbances for all L' . For < a' < a' c , 
where a' c ~ 5.7, all traveling waves are unstable. For a' > a' c , a stable region appears (see 
Fig.[S^e)). Subsequently, for any value of a' > a' c the stable region becomes more apparent. 

This result can be compared with the numerical stability results shown in Fig. Ufa), 
where a = 1/2. The stability condition a' > a' c ~ 5.7 (where a 1 = ct/y/r) corresponds to 
r < (a/5.7) 2 = 0.0077, showing remarkably good agreement with the upper limit of the 
stable region in Fig. Efa). 

If instead we consider the special case a' = 0, with 0' > 0, we find a broadly similar 
picture, in that all traveling waves are unstable when /3' is small, but some eventually 
stabilize, once (3' is sufficiently large. From Fig. [9] it is apparent that the two stability 
boundaries R and T intersect, coalesce, then lift off from the q' axis as /?' is increased. 
Ultimately they re-attach to the q' axis, when /3' = (3' c , where f3' c ~ 5.06 as shown in 
Fig. [9fg). For 0' > (3' c , there is a region of stable traveling waves. 

Let us now express the results above in a form more illuminating for comparison with our 
earlier numerical secondary stability calculations (Sec. IHIj) . As an example, we set a = 1 
and j3 — and consider the limit of small L, looking for regions of stable waves as r 2 is 
varied. 
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FIG. 9: The stability boundaries for traveling waves according to the amplitude equations (|24[)— 
(1261) . in the special case a' = 0, for different values for p: (a) /?' = 0.2, (b) /?' = 0.476, (c) ft = 0.6, 
(d) p = 2, (e) 0' = 4.5, (f) = 5, (g) f3' = 5.0607, (h) f3' = 6. See the text for more details. 

From rolls are unstable as long as q > llr 2 /12 + {) 2 /1152. If q < llr 2 /12 + £> 2 /1152, 
then 0i is purely imaginary and hence 02 must be considered. From fl29l) we have definite 
instability if r 2 > 144g/91. In addition to these rather blunt conditions, the sign of 02 
must also be considered in order to determine the stable region. Figure [10] shows the curves 
q = llr 2 /12 + -0 2 /1152 (solid line), q = 91r 2 /144 (dashed line) and 2 = (dotted lines). 
Any region of stability must lie between the solid and dashed lines. After checking carefully 
the signs of the eigenvalues, we find that the stable region (indicated by the asterisks in 
the figure) lies between the two dotted lines in the upper and lower parts of the graph, and 
between the dotted and solid lines for a small range of intermediate values of r 2 (see Fig. [TO]) . 
Although they appear almost parallel in Fig. HOT a), for large r 2 , as in Fig. [TOT b). the two 
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FIG. 10: The region of secondary stability of traveling waves is marked with asterisks; for details 
refer to the text. The solid line shows where a\ is purely imaginary; to the right of this line, the 
traveling waves are certainly unstable. The dashed line shows where q = 91r 2 /144; to the left of 
this line, traveling waves are also certainly unstable. The dotted lines show where cr 2 = 0. 

sides of the secondary stability region are no longer approximately parallel. 

Note the qualitative similarity between the shapes of the stable regions in Fig. [10] and 



The question remains of whether or not this stable region extends to indefinitely large 
values of r 2 . To investigate the large-r 2 behavior of the stability region, we consider large r 2 
with q = 0(r 2 ), motivated by the observation, from Fig. [TOT b). that stable rolls lie in some 
region between straight lines in (q,r 2 ) parameter space. In this limit the stability condition 
from ( 1281) simplifies to q < llr 2 /12, while cr 2 = 91/2 — 72g/r 2 + 0{r 2 ). Hence we can 
conclude that the region of stable waves for small L and large r 2 is 



In summary, the results of this section show that when a and /3 are O(e), there can be a 
narrow region of stable traveling waves near k — 1, and that there is no upper limit on the 
size of r 2 allowing stable rolls. 

For even smaller values of a and (3, of order e 2 , we have checked that a and f3 do not 
appear in the leading order amplitude equations, so in that case all traveling waves are 
unstable, as in the non-dispersive case. 



Fig. El 



91r 2 /144 < q < llr 2 /12. 



(30) 
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B. Traveling waves close to the marginal curve 



We now turn to the second case in which a^q may be small: in the region close to the 
marginal stability curve. Following an analysis similar to that for the dissipative Nikolaevskiy 
equation |4|, we find that, in contrast to the dissipative case (in which a narrow region of 
stable rolls exists close to the marginal curve Q), here all traveling waves are unstable near 
the marginal curve. 



VII. NUMERICAL SIMULATIONS OF THE DISPERSIVE NIKOLAEVSKIY 
EQUATION 

To illustrate some of the consequences of the results of the preceding sections, we have 
carried out numerical simulations of the dispersive Nikolaevskiy equation, using a pseu- 
dospectral method, with exponential time stepping of which a small sample are pre- 
sented here. The initial condition is taken to be a traveling wave with a given wavenumber 
k (approximated as a cosine of the amplitude given by ©), plus small random noise, and 
the domain size is D = lOOir/k. Figure ITT1 illustrates in order strong (a), intermediate (b) 
and weak dispersion (c)-(e). The values of a, (3 and k are chosen in each case to correspond 
to traveling waves that are predicted to be unstable by the asymptotic analysis. 

Figured]] (a) shows the case of strong dispersion, with a = 2 and (3 = 1, and wave number 
k = 1 and r = 0.01. The stability analysis of ( ITT)) predicts that the rolls are unstable, since 
a = 2 and = 1 lie in the unstable region in Fig. |5j the numerical simulation agrees with 
this asymptotic result. 

In Fig. [TT] (b) an example of intermediate dispersion is simulated, where a = 2e 3//4 and 
(3 = e 3 / 4 , r = e 2 and k = 1 + eq, for q = 0.2 and e = 0.1. It is known from the asymptotic 
results of Sec. [V] that a and (3 being 0(e 3 / 4 ) with wave number k = 1 + eq will result in 
unstable traveling wave solutions, which agrees with the simulation shown in Fig. [TT] (b). 

To show the effects of weak dispersion with wave number k = 1 + e 2 q we take r = 0.01 
and e = 0.1. Figure [TTJ (c) shows the case a = 2e, (3 = 0, q = 0.87, while in Fig. [TTJ (d) 
the parameter values are a = 0, (3 = 5e, q = 2.5. Rolls should in each case be unstable, 
according to the analysis of Sec. IVI A) and this is confirmed by the numerical simulations. 

Figure [TT] (e) represents weak dispersion, with a = e and (3 = 0. The wavenumber is 
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FIG. 11: Snapshots of the numerical solutions of ([2]). Parameter values are: (a) a = 2, f3 = 1, 
r = 0.01 and k = 1; (b) a = 0.3557, (3 = 0.1778, r = 0.01 and k = 1.02; (c) a = 0.2, /3 = 0, r = 0.01 
and k = 1.0087; (d) a = 0, = 0.5, r = 0.01 and k = 1.025; (e) a = 0.25, /3 = 0, r = 0.0025 and 
& = 1.00125. The times of the snapshots are indicated in the insets. 

k = 1 + e 2 g and r = e 2 r2, for e = 0.25, q = 0.02 with r2 = 0.04. These values of r2 and 
q lie in the unstable region given in Fig. [TO], and the simulations support this prediction of 
instability. 

VIII. CONCLUSIONS 

We have examined the stability of spatially periodic solutions to the dispersive Niko- 
laevskiy equation, which is the original model introduced by Nikolaevskiy [lj for seismic 
waves. The reincorporation of dispersive effects stands in contrast to most studies subse- 
quent to Nikolaevskiy's paper. We have shown how the instability of all spatially periodic 
solutions at the onset of pattern formation in the more-often treated, nondispersive version 
is modified by the presence of dispersive terms. Our results have been achieved through both 
a numerical calculation of the secondary stability boundary for the traveling wave solutions 
and an asymptotic treatment of three particular scalings in e for the dispersive terms. The 
secondary stability diagrams ( "Busse balloons" ) can be rather complicated, and can depend 
sensitively on the size of the dispersive terms. 

Our consideration of the case a, (3 = 0(1) can be interpreted as giving information about 
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the bottom of the secondary stability diagram obtained in (k, r) parameter space for fixed 
a and 0. Two cases were found: either all traveling waves are unstable at the bottom of 
the diagram, or there is a symmetrical, Eckhaus-like region of stable traveling waves, right 
down to onset at r = (although the width of the region of stable rolls does not stand in 
the usual Eckhaus ratio to the width of the existence region of rolls). 

The separate analysis for smaller values of a, (3 can be interpreted as shedding light on 
the upper parts of the fixed-a,/3 stability diagram in k, r parameter space. We have shown 
that for small a, (3, a narrow region of stable waves may exist near k — 1. However, beyond 
the range of validity of the asymptotic analysis, the numerical stability results show the 
complicated nature of the secondary stability boundaries, so we are unable to draw any 
significant general conclusions about the form of the secondary stability diagram, limiting 
ourselves to some specific examples. Things are further complicated by the fact that rolls 
predicted to be stable by the asymptotics may in fact turn out to be unstable when the full 
numerical calculation is performed, since the asymptotics concerns only long-wavelength 
instabilities, and other, short-wavelength instabilities may turn out to be present. 

In this paper, we have said little about the behavior of time-dependent solutions of the 
dispersive Nikolaevskiy equation. However, it appears from our numerical simulations that 
when all waves are unstable, chaotic states are found that have a similar behavior to that 
found in the non-dispersive Nikolaevskiy equation {5- 7]. 
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